The Adjoint Problem in the Presence of a Deformed Surface: the Example of the 

Rosensweig Instability on Magnetic Fluids 
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The Rosensweig instability is the phenomenon that above a certain threshold of a vertical mag- 
netic field peaks appear on the free surface of a horizontal layer of magnetic fluid. In contrast 
to almost all classical hydrodynamical systems, the nonlinearities of the Rosensweig instability are 
entirely triggered by the properties of a deformed and a priori unknown surface. The resulting 
problems in defining an adjoint operator for such nonlinearities are illustrated. The implications 
concerning amplitude equations for pattern forming systems with a deformed surface are discussed. 

PACS number(s): 02.30.Tb, 47.20.-k, 75.50.-y 



I. INTRODUCTION 



Magnetic fluids (MFs) are stable colloidal suspensions of ferromagnetic nanoparticles dispersed in a non-magnetic 
carrier liquid. The nanoparticles are coated with a layer of chemically adsorbed surfactants to avoid agglomeration. 
The behaviour of MFs is characterized by the intricate interaction of their hydrodynamic and magnetic properties 
with external forces. This complex interaction causes many fascinating phenomena, as the 'negative viscosity' effect 
and the Weissenberg effect (for a review see jjj) or as the labyrinthine instability and the Rosensweig instability ||]. 
The latter instability occurs when a horizontal layer of MF with a free surface is subjected to a uniform and vertically 
oriented magnetic field. Above a certain threshold of the magnetic field that surface becomes unstable, giving rise to 
a hexagonal pattern of peaks ^jj^]. Despite the fact that the Rosensweig instability has been known for many decades, 
some aspects have been addressed only recently: the inclusion of the fluid viscosity [Q and of a finite layer thickness 
[^J, the hexagon-square transition |(| or the wave number selection problem ^,^]. 

The quantitative comparison of theoretical and experimental results is presently limited to linear aspects of the 
Rosensweig instability. Convincing quantitative agreement is found for the wave number of maximal growth H , for 
the parametric stabilization of the Rosensweig instability || , and for the wave resistance in magnetic fluids For 
the nonlinear aspects, comparisons are restricted to the detection of the same qualitative features. Stability regions 
for hexagons and squares are found in the theory [^0[-|l^] and in the experiment ||. The experimentally observed 
increase in the wavelength of the emerging hexagonal pattern ft] is detectable as well in a theoretical analysis (l6j . 
The main reason for lacking quantitative comparisons is the difference between the susceptibility of the experimentally 
used fluids and the susceptibility up to which the nonlinear analyses are valid. Finite amplitudes for all three patterns 
of hexagons, squares, and rolls are given for susceptibilities x smaller then 0.41 for an infinitely thick layer (see Figs. 3 
and 6 in [fl6||). But the experiments were performed with a magnetic fluid of Xcxp = 1-4 ((J. Pattern selection studies 
.with fluids of different susceptibilities have not yet been carried out. 

Weakly nonlinear analyses of the Rosensweig instability were done by means of an energy minimization principle 
""IT by methods of functional analysis p^l5[|, by a generalized Swift-Hohenberg equation |l7|,[l8), and by a 
tiple-scale analysis ji~9|,^0) . The first two approaches are suited only for static problems and the Swift-Hohenberg 
J> | equation lacks coefficients containing the fluid properties and the geometry of the system. Consequently, amplitude 
k>( ' ec l ua ti ons stemming from a multiple-scale analysis based on the fundamental hydrodynamic equations are customized 
rS | for the study of static and dynamical problems as well as for the quantitative comparison with the experiment. The 
^ .standard route of the multiple-scale analysis is modified in |l^,|2(| to circumvent Fredholms theorem, i.e. the definition 
of an adjoint operator. In the present paper a multiple-scale analysis is presented which involves the expansion of all 
physical quantities at the deformed surface. The resulting problem in defining an adjoint operator and the subsequent 
consequences are the main purpose of this paper. It is organized as follows: the system and the relevant equations 
of the problem are displayed in the next section. Based on the governing equations and th e bo undary conditions the 



different character of the nonlinearities of the Rosensweig instability is emphasised (Sec. II C). Whereas the linear 
problem is shortly recapitulated in Sec. Ill, the adjoint problem is addressed in detail in Sec. [V. In the final section 



the problems concerning amplitude equations for pattern forming systems with a deformed surface are discussed as 
well as open questions and further prospects are outlined. 
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II. SYSTEM AND EQUATIONS OF THE PROBLEM 



A horizontally unbounded layer of an incompressible, nonconducting, and viscous magnetic fluid of finite thickness d 
and constant density p is considered. The MF is bounded from below [z — —d) by the bottom of a container made of a 
magnetically impermeable material and has a free surface described by z s = y, t) with air above. The electrically 
insulating fluid justifies the stationary form of the Maxwell equations, which reduce to the Laplace equation for the 
magnetic potentials «&W in each of the three different regions. Upper indices denote the considered media: (1) air, 
(2) magnetic fluid, and (3) container (Fig. [I]). 
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FIG. 1. Sketch of the system. 



It is assumed that the magnetization M' 2 ' of the magnetic fluid depends linearly on the magnetic field H^ 2 ', 
M< 2 ) = (p r - 1)H( 2 \ where p r is the relative permeability of the fluid. Additionally, the magnetization is considered 
to be a linear function of the density p^ which results in the usual form of the Kelvin force Fk = p>o (M( 2 )grad)H( 2 ) 
0. The system is governed by the equation of continuity, the Navier-Stokes equations for the magnetic fluid, 

divv (2) =0, (1) 
p^d t v^ + p ( 2 ) ( v ( 2 )grad)v( 2 ) = -gradp< 2 > + // 2 > Av^ + ^(M^grad)^ 2 ' + p^g , (2) 

and the Laplace equation in each medium, 

A$W = . (3) 

The scalar magnetic potentials are defined by — ~ grad'&W. The velocity field in the MF is denoted by v' 2 ), the 
dynamic viscosity by p( 2 \ the pressure field by p^ 2 \ and the acceleration due to gravity by g = (0, 0, —g). The first 
three terms on the right-hand side of Eq. (||) result from div T ^ , where the components of the stress tensor T ^ 
read @ 

Tf = j- P (2) - ^o^} a« + nf^Bf + p^ ( diV ? + 9j 4 2) ) ■ (4) 

M< 2 \ H {2 \ and B (2) denote the absolute value of the magnetization, the magnetic field, and the induction > in 
the MF. Since the final arrangement of peaks is a static configuration, the velocity field in the MF is set to zero, 
v' 2 ) = 0, and the surface depends only on the horizontal spatial coordinates, z s = ((x, y). Applying all assumptions, 
the static form of the Navier-Stokes equation is 

= -gradp( 2 > + f (A* - 1) grad(tf < 2 >) 2 - p^ge z , (5) 

where e 2 = (0,0,1) is the unit vector in z-direction. The governing equations have to be supplemented by the 
appropriate boundary conditions. These are the continuity of the normal (tangential) component of the magnetic 
induction (magnetic field) at the top and bottom interface, 

n (2,i) . ( B (i) _ B ( 2 )) = 0, n^ 1 ) x (H« - H( 2 )) = atz = z s , (6) 

n (3,2) . / B (2) _ B (3)j = Qj n (3,2) x / H (2) _ H ( 3 ))= at2 = -d, (7) 
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and the continuity of the normal component of the stress tensor across the free surface 

/<y>(i)_<y>( 2 )\ ( 2 ,i)_ (TKn (2,i) =0 atz = z s . (8) 



In (||) the surface tension between the magnetic fluid and air is denoted by a, the curvature of the surface by 
K = divr/ 2,1 ), and the unit vector normal to the MF surface by 

n (2,i) _ grad [z - ((x, y)} = {-d x (j, -d y £, 1) (Q) 
~ |grad[z-C(z,y)]| V 1 + (d*0 2 + (d y C) 2 ' 

The upper index (2, 1) at the unit vector indicates that n/ 2,1 ) points from medium 2 towards medium 1; analogous 
for the normal vector n^ 3 ' 2 ) = (0, 0, 1) (see Fig. [l]). The difference of the tangential components of the stress tensor is 
identically zero because of the continuity of the magnetic fields and inductions in (J6|) . Neglecting the influence of the 
air pressure with respect to the fluid pressure, p^ ~ 0, and using = 0, one finally gets from Eq. (||) 



Mo 
2 



M (2) n (24) j 2 + p {2) - oK = atz = z s , (10) 



where p^ is the solution of Eq. (Eh. 



A. Basic state 



As long as the applied spatially homogeneous magnetic field perpendicular to the surface, 

H« =Hg } = (0,0,i4 l) ) BW = b£9 = (0,0,B$) = Mg } = (0,0,M^ 2) ) , (11) 

is below a certain strength, the system is in its basic or ground state. This state is given by the plane surface z s = 0. 
The corresponding solution for the fluid pressure is 

J»g>=-p<V-f (A4 2) ) 2 , (12) 

where the constant resulting from (j|) was determined by ( |To) ) and Hq' — Hq^\ was used. Here | denotes the 

evaluation a1 
respectively. 



evaluation at the plane interface z s = 0. The boundary conditions (^, 0) are fulfilled by Bq^ = Bq^ and Bq = Bq , 



B. Small disturbances and their expansion 

In order to study the stability of the basic state and the pattern selection problem in the weakly nonlinear regime, 
small deviations from the basic state are considered 

z s = + C, H« = H$ + h« , B« = B$ + b« , 

M( 2 ) = Mg } + m( 2 ) , *W=$W+0(O i p W =p^+nW . (13) 

Dimensionless quantities are introduced, where physical quantities are denoted by a hat in the rest of the paper. All 
quantities associated with the magnetic field are scaled by the critical magnetic induction in the limit of an infinite 
thickness of the MF layer, B c .oo- This preference is based on the fact that magnetic inductions can be directly 
measured by Hall probes. Furthermore, it is assumed that the deviations are proportional to the applied external 
induction. Thus the following scaling is used 

Bg ) - B cxt B c ,oo , b w = b«B cxt B c ,oo , (14) 



M« = Bcxt ^-il4,oo , m« = m W Bext (^_i) Bci00 , (is) 

AoMr AoMr 

Hg> = B cxt ^ , h« = h«i? oxt ^_ , (16) 

AoA*r MoMr 

^(2)^,(2)^2 i = J-. (17) 

Mr (Mr - 1) fcc.oo 
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For an infinite thickness of the layer, the critical induction and critical wave number, respectively, are 

% 



b: 



2/t fJriPr + l)Vp (2) 0-g 



(2), 



(18) 



The dimensionless quantity -B C xt in B ex t = (0, 0, -B C xt) measures the strength of the applied external induction in units 
of B coo . 

With this scaling the solution of the Navier-Stokes equation (JsJ) becomes 



(2) V 

1-ry 



2h^ + (h^) 



with a yet unknown constant c. The Laplace equation (|^) in each medium is 

A0« = , 

and the boundary conditions (|[ (?], |l^) arc 

&tz = -d = -<t> {2) + (3) , 



at z = — rf 
at z = £ 
at z = £ 



1 — 7J 

2r? 



c + 



(1) - ^(2) 



V 1 - v 

■d.4P> ' 1 + 11 



(2) 



1-ry 



9, 



at z = C = - 



+ 



2^ [(^cr+iMD 

(! + »?) [1 + {d x Q 2 + (dyC) 2 } 
2B e 2 xt (l-77) / 2 V h^nP' 1 ) 



h( 2 )< 



+ 



?7(1 + rf) 



2r? (h( 2 ) n ( 2 ' 1 )] 
(l-»7) 



(h (2) )^ 



r, 2 (l + ry) 



The curvature of the surface is 



K = 



$xx C ^yy C 



(d X () 2 d XX C + 20 X ( <)„>■ 8 X yC + (OyO 2 8 yy C 



y/1 + (d x C) 2 + (dyC) 2 

and the widely used quantity 



[1 + {d x Q 2 + (dyCY 



i 3/2 



/i r + 1 



(19) 

(20) 

(21) 
(22) 

(23) 
(24) 



(25) 



(26) 



(27) 



was introduced (7 = (3/4)ry in Q, ©0). 

Since the disturbances of the magnetic field are located in the vicinity of the magnetic fluid, they have to disappear 
as z goes to infinity, i.e. the magnetic field disturbances have to fulfill two more conditions: 



at z = -00 = (3) 
at z = 00 = c/> (1) 



(28) 
(29) 



In summary, the system is gov erned by three differential equations of second order ( |20| ) and the corresponding six 
boundary conditions (pl|-p4|, [28, 29 ) for the magnetic potentials. The remaining equation p^) is a differential equation 
for C(x,y) whose solution determines the shape of the surface. 
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To derive solutions with a finite but small amplitude, the Eqs. (|20|-|2, r 
control parameter B ex ± and the various physical quantities are expanded 



B. 



are solved perturbatively. The external 

(30) 
(31) 

where X can be any of the quantities ) £ ; and c. e is an expansion parameter with respect to the magnitude of 
the surface deformation (0 < e <C 1). The boundary conditions ( p3f]25| ) have to be evaluated at an a priori unknown 
surface £(x,y). Therefore it is convenient to expand all potentials p ob the surface in terms of cf>^ at z = 0, 



cxt 

X 



B z c +eBi+e z B 2 + 
eX +e 2 X 1 +e 3 X 2 



0« (x, y,z = = <P W (x, y, 0) + d z <f>w (x, y, z) 



- n - ,/,(*) 



(i)t 



■C + -d zz ^{x,y,z) 





c 2 + 



(32) 



and similarly for all space derivatives of 4> . The chosen expansion of the external parameter follows the way as it 
is known from thermal convection problems of ordinary fluids [ p2|]23j | or magnetic fluids pi| and from the Benard- 
Marangoni problem p5|-|2^]. The expansion^ 30, |l|) describes a situation slightly above the threshold of the linear 
instability and differs from those used in [[19|2^], where the magnetic field H, acting as external parameter, is not 
expanded. The consequences of this difference will be discussed in Sec. [v| 



C. Different nonlinearity 

Before going into the details of the expansion procedure, it is illuminating to analyze the set of Eqs. (^0|-p5|). The 
differential equations ( poj ) and the boundary conditions at the bottom of the container (^, p2| ) are linear equations. 
Only the three remaining boundary equations (p3|-p5|) contain nonlinear terms. In contrast to almost all classical 
hydrodynamic systems, the nonlinearity stems from the boundary conditions and not from the differential equations 
which describe phenomena in the bulk of the system. For example in convection systems the nonlinear contributions 
are caused by the advective term, i.e. by a bulk term, in the Navier-Stokes equation and in the equation of head 
conduction. For the Rosensweig instability the nonlincaritics are caused by the deformed and unknown surface £ at 
which the potentials have to be calculated. Thus each linear term in tp^ and 4>^ in Eqs. ( p3f]25| ) generates nonlinear 
terms via Eq. (^). A further nonlinear contribution comes from the Kelvin force density in which a term quadratic 
in the potential (f>^ appears, to determine again at the deformed and unknown surface. This qualitatively different 
character of the Rosensweig instability makes it fascinating as well as hard to treat the instability with the classical 
methods. One indicator of the essential role of the boundaries is the fact that the external parameter B ext appears 
only in the boundary condition (p5|). 

Performing the expansion according to (pb]-|32]) in (|2(] ||^) generates a hierarchy of linear equations for Uo, U\, and 

U 2 



L Uo 
L U 2 



0, 



-I^Uo+N^UoM) , 
-L 1 i7 1 -L 2 C/o + A r (C/o,C/i) 



NiU^Uo) + N(U Q , Uo, Uo) + N 2 (U , U ) 



by matching the three lowest powers of e. The vector 

U = , <p&, 0« 0« | Q , (2) | , C) T = eUo + s 2 U t + e*U 2 + 



(33) 
(34) 
(35) 

(36) 



is the augmented state vector, Lo, L\, and L 2 are linear operators, and the vectors N(-), Nx(-), and N 2 (-) contain 
the nonlinear contributions. With the above mentioned, it is clear that in each order of e different nonhomogeneous 
boundary problems have to be solved. This again indicates the different quality of the problem considered here 
in comparison to classical convection problems. For the boundary conditions, terms involving (f>^, 4>^ 2 \ and their 
derivatives have to be evaluated at z = wherefore the vector U contains the potentials tfy^ and tft^ at z = 
additionally to the potentials itself. Similar state vector were used for the Benard-Marangoni system or for 

the electroconvection problem in a thin film |28[ . 



III. LINEAR PROBLEM 



The linear stability problem and the corresponding adjoint problem play a fundamental role in the nonlinear analysis. 
Since the linear problem has been considered in details elsewhere [Q,[|,[7 29 3lJ only some of its more important aspects 
are recapitulated here. 
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Inserting the ansatz (B0[]32) into (|20H23) gives at order 0(e) 



= A0« , 



at z = — d 
at z = —d 

at z = 

at z = 

at z = 



= -. 



o = -I±^ 2 > 

1 -ry 



1-7, 
(1) 



= -d z <f 

: 



l+?7 



1-7, 

2(1 - V )B? 



,.. ; , L ^4 X) ~ CO + (3„ + dyy) CO + (l 2n l 2 ^ C0 
77(1+7?) 7,2(1+77) 



The equations ( |37T|ll| ) are solved by 



v ( 1 + Tl)e k(d + z) ^ 



Coe 



(1) 7,(1+7,) [ e Hd-z) _ ve -k(d +Z ^ ^ 

----- — — <,oe 



_ ?7 ) [ e fcd _ ^g-fcd] 

Co - Coe lkr , 



(37) 
(38) 

(39) 
(40) 
(41) 
(42) 

(43) 

(44) 

(45) 
(46) 



with k = (k x , k y ), k = y/k% + fcjj, r = (x,y), and Co the amplitude of the surface deformation. Inserting this solution 

into ( |42|) leads to the known dispersion relation for an inviscid magnetic fluid of finite thickness in dimensionless 
quantities [ff^Mpf} 



1 + fc 2 - 2B 2 r k 



Jid 



-kd 



^kd 



rj 2 e kd 



= 0, 



(47) 



providing that the unknown constant Cq is set to zero. By solving equation (fl7j ) numerically for a given thickness d, 
one gets the height dependent critical values B c (d) and k c {d) at the onset of the instability. 

The linear operator Lq together with the corresponding augmented state vector Uo of Eq. (|3^) are defined on the 
basis of Eqs. (M HjfP): 



(48) 



L U = 
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d xx 
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IV. ADJOINT PROBLEM 



To proceed with the derivation towards an amplitude equation, one has to solve the linear inhomogeneous equations 
for Ui and U2, Eqs. (|4, 35). According to Fredholms theorem |Q, they have a solution if and only if the right hand 
side is orthogonal to the zero space of the linear operator L . If Uq is an eigenvector of the adjoint linear operator Lq 
with zero eigenvalue, then multiplying Eq. (Q) from left with Uo reads (analogous for Eq. (|35|)) 



< U , L x Uo > + <U Q , N^Uo, Uo) >=< U , LqU, >=< LjE7 , U Y >= 



(49) 



The scalar product is denoted by < , >, where the explicit form depends on the considered problem. Applying these 
solvability conditions, by which B\ and Bi (see Eq. (|30|)) can be expressed, leads to the general amplitude equation 
for a pattern selection problem. 
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In Sec. [I C the different character is described of the pure 'surface-nonlinearities' occuring in the Rosensweig 
instability. The question arises whether the established methods for 'bulk-nonlinearities', as for example for the 



Rayleigh-Benard convection, will work here as well. Motivated by 
and the plain surface appear, the following scalar product 



where nonlinearities triggered by the bulk 



<U ,U > = lim 



A/2 


A/2 




/ dx 


/ dy 


[/: 


1-1/2 J 


-1/2 





dz 



. dz$ ) *<j>^ + 
1 — V J-d 



dz 



+du* 0A | o (1) 1 + e< 5 1 Q <t>& | Q + fu* Qfi | o C 



(50) 



is used. Uq is chosen since Eqs. (|34|, |35|) have to be multiplied by Uq. The first three terms in the square brackets 
in Eq. ([SO]) are contributions from the volume, whereas the last three terms are contributions from the surface. The 
possible prefactors d, e, and / schwach nichtlinear as well as the surface components of the adjoint vector Uq, uo,4, 
mo. 5, and Mo. 6, have yet to be determined. From [^6|,^] it is known that the surface components of the adjoint state 
vector are not just simply the adjoint components of the state vector. Using the identity 



$>A<t>M = dx d X ^ - <£«1 + dy dy^ - 



d 



XX'f'Q 



+ d yy4>o 



00* Ai) 



a. 



(51) 



and the conditions ( |2l| , p2| , £ 

< U ,L U > 



one has after partial integration 



lim i 



1/2 



dx 



1/2 



1/2 



1/2 



dy 



dzA^*^ 



1 



n 



+ / < teA^>V 1 J-^|_ 00 ^W|_ 00 +(^*-^*)| 



1 — t] 

1(3)* 



dz A0 O ' 



(2), 



-a. 



>ZH>Q 



1 



n 



</> (3) l 



-d 



(52) 



From Eq. 



it is concluded that the linear adjoint system is governed by three differential equations of second order 



for the adjoint potentials 



According to Eq. (152) the four conditions to be fulfilled by these adjoint potentials are 



at z = — oo 
at z = —d 
at z = —d 



= <A 



(3) 
'0 



- 1(2) 



= - 



1 



V 



1-7? 



;( 2 ) 



zV 



at z — oo I 

since only terms on the surface should remain. The residual terms 



(53) 
(54) 

(55) 
(56) 



R = d z ^ 2 
+d z d> 



1 



2r? 



1 



1 



du. 



1 + T] 

T 



■ eu 



0.5 



1 



l 



n 



2V0 



dUn 



2(1-7?) a 



0,4 



(0* 



9 



7?(l + 7?)' 
1) /«$,. 



^ u 0,4 



(57) 



are bounded to contain only the two remaining boundary conditions and one further term, all to be determined at 
z = 0, for a proper definition of Lq. Comparing the terms in R and the surface components of U (sec Eq. (jffij)) . 
one realizes that the first and third term in Eq. (|57|) have to vanish. As a consequence of the immediate choice 



4 2 \ 



the sixth component of the adjoint vector Uq, uq q, is a linear combination of 0q and 



4 2) , 



i.e. a linear 



combination of (j) and tto,5- Certainly a result contrary to the condition that the set of variables in Uq have to be 
linearly independent. Testing the choice 

du A = fc<W, e "o,5 = 



»o 2) 



-?? 2 (l + r?) 
(1 - rtYBl 



Co 



(58) 
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results in 



< U ,L U > 



lim — 

L^oo I 2 



A/2 


A/2 




/ dx 


1 dy 


[/-' 


1-1/2 J 


-1/2 





dz a4 1} * 



v 2 (i + v) 

(1 - rtfBl 
(1 - J?) 



^-v) B c 5 7(1), 



Co* + 4^* 



(21* 



- V 



J yy 



V J- 
o 

- l)Co 



I) 



dzA^*^ 



(59) 



where <9, 



6q I q was use d- This identity is justified because of the Laplace equation for (f>jp (see last 



'zVO |o 

^-integral in Eq. (|52|)) and the condition (|56|). After inserting ( p*8| ) into Eq. (f50|), it becomes clear that Eq. ( |59| ) has 
almost the proper form in order to define the adjoint operator Lj via < Uo, LqU >=< LqUq, U >. The last remaining 
but unsolved problem is posed by the third surface term in Eq. (|59|). Instead of fc^^jo, the expression — d z (f)^\o 



appears. Both terms are only equally in the first order of the expansion, where cj^p ~ e~ kz (see Eq. (|45[)) . In higher 

orders of the expansion, the functions <j>^ and <fi^ containing wave vectors which are linear combinations of two and 

three, respectively, wave vectors of the basic modes. As a result z-dependences will appear as cj£p ~ e _ l k ™ ±km l 2 and 

4>^ ~ e _ l k '* ±km±ki l z , respectively. The absolute value of the resulting wave vector is usually not equal to k. Thus 
the task remains that for a nonlinear analysis a linear adjoint operator should be definable which is valid for a wider 
set of function as e~ kz . 

The generic problem of a nonzero, deformed surface is illustrated by calculating the scalar product of the augmented 
state U vector and its adjoint one U 



<U,U> = e 2 < U , U > +e 3 (< U , U x > + < U u U >) + G>(e 4 ) , 



lim -rx 



A/2 


A/2 




/ dx 


/ dy 


[/- 


1-1/2 J 


-1/2 





dzft 3 >cf>® 



1 + 77 



V J-d 



dzft 2 >^ + / dzftV*^ 



(60) 



Since U and U contain the complete information of the disturbed state, the unknown and deformed surface £ appears as 
a bound at two integrals and at the surface contributions. Whereas the expansion of the three surface contributions 
can be performed accordingly to Eq. ([32]), the expansion of the two integrals involving the deformed surface £ is 
accomplished as 



C dzft 2 >^=e 2 f dz<p (2) *<}> {2) +e 3 

-d J-d 



Co0o IcA lo + 



dz 



— d 



2) 



+ 0(e 4 ) 



(61) 



The analogous expansion can be applied to the second integral involving £. Matching the two lowest powers of e in 
Eq. (|60|), one gets 



< U ,U > = lim — 



A/2 


A/2 




/ dx 


1 dy 


if. 


1-1/2 J 


-1/2 





dz 



WW 



+du*o A \ 4> { o 1] \o + eu h\o ( t )i o ) \o + / M o,eloCo 



V 



1 — 77 



dz 



WW 
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< U Q , Ux > + < Ux, U Q > = lim \r 
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(63) 



As discussed for Eq. (|59|), the problems are appearing beyond the first order of expansion. The form for the scalar 
product of < Uq, Uq > is the expected one (Eq. (|62|)). In the second order expansion, the left hand side of Eq. ( |63| ) 
lets await only products of functions which belong to the first and second order of expansion. But the expansion 
on the deformed surface generates additional terms, where each one is formed by three functions of the first order 
expansion. These additional terms can neither be assigned to < Uq,Ux > nor to < Ui,Uq > by simple arguments. 
Similar dilemmas are present for the third and any higher expansion order with many more not-assignable terms. 
Neglecting surface deformations is certainly not an option for the solution of this problem. 

It becomes now apparent why the derivation of an amplitude equation in 26 - 2£ 
surface succeeded. 



with the approximation of a plane 
The first order expansion on a deformed surface is equivalent to that approximation of a plane 
surface. Therefore the problems of not-assignable terms stemming from higher orders in the expansion do not occur. 



V. DISCUSSION 



The detailed analysis of the last section showed that the definition of a linear adjoint operator for pure surface- 
nonlinearities has not yet been accomplished. Two attempts which led to unsatisfactory consequences were described 
here. They point to the fact that in the known literature no derivation can be traced in which an expansion on 
the deformed surface was taken into account and a linear adjoint operator was defined. The latter is necessary to 
determine higher order terms of the external parameter; for the system considered here these are Bx and B2, see 
Eq. (32). This fact is surprisingly since pattern forming systems with a deformed surface belong to the set of classical 
hydrodynamical systems. 

For the nonlinear analysis of instabilities in MF, a very similar expansion route was followed in |il|^0],[33| . Re- 
markably, only very shortly the existence of an adjoint problem is mentioned in |^,^]. No details with respect to 
the definition of a scalar product and a linear adjoint operator were published. The same expansion as in Eq. ( ]3l| ) is 
used, but the external parameter, the magnetic field H, is not expanded accordingly to Eq. (^0|). The latter step is 
a inconsistency since the multiple-scale expansion comes along with the expansion of the physical quantities and the 
external driving parameters as the Rayleigh number [p2|-p4| or the Marangoni number J25|-^7j . The gain of the lacking 
expansion is that higher order terms of the external parameter have not to be determined. The solvability condition 
in second order postulated in fl9| , |20|| demands that the amplitude has to have a nonvanishing derivative everywhere 
with respect to the slow spatial variable. That means that the amplitude is a strictly monotonously increasing or 
decreasing function which is a very special type of solution. Selection problems between regular patterns can not be 
tackled with such a type of solution. 

The problems caused by a deformed surface for the derivation of an amplitude equation are not unique to instabilities 
of magnetic fluids. The analysis of the Benard-Marangoni convection or of parametric surface (Faraday) waves is 
confronted with similar difficulties. For the Benard-Marangoni convection the surface deformation is often disregarded 
in fluid-gas systems [ ^6|j3^ 1 and two-fluid systems [^tJ. In j3^J|(| a surface deflection was imitated by a nonzero 
Crispation number, thus avoiding the expansion of any quantity on a deformed surface. In pq| the expansion of 
the layer thickness was tuned in such a way that eigenfunctions could be used which correspond to the case of an 
undeformed interface. 

In the nonlinear analysis of Faraday waves an explicit evaluation of quantities on the free surface is rare. In |37|| 
the driving parameter is not expanded, similar to ]l9| , p0| , ^3[ . No expansion on the surface is performed, but two 
different amplitudes are introduced for the solution of first and third order. The adjoint operator is circumvented in 
this way. The two different amplitudes are just the two quantities needed to attach use to the solvability conditions. 
Under the constrains of an inviscid fluid and an irrotational flow field, Milner [B8J derived an amplitude equation for 
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parametrically driven capillary waves. An expansion of the velocity potential on the surface was included, but an 
adjoint problem was not formulated. 

The overall picture is that several circumventing routes were used in order to avoid the explicit evaluation of 
quantities on a deformed surface and the formulation of a linear adjoint operator. The used alternatives were hardly 
motivated which is why the real reasons behind the search for alternatives are not publicly known. Instead of using the 
solvability conditions to determine higher order terms of the external parameters, they were converted to constraining 
conditions for the amplitudes of the expanded physical quantities. This presents a rather unsatisfactory situation. 

The deeper reasons for this unsatisfactory situation are disclosed in this article. The adjoint problem is presented 
in detail for nonlinearities purely triggered by a deformed surface as in the case of the Rosensweig instability. By 
expanding the physical quantities and the external parameter, the knowledge of the linear adjoint operator is essentially 
to proceed towards the nonlinear amplitude equation. The main and still unsolved problem is the proper definition 
of the linear adjoint operator. Since this problem mounts a principal barrier and is not mentioned as a reason for 
the search of alternative approaches, it is presented here even if no solution can presently be offered. This situation 
entails one foremost question: Is it possible to apply the concept of Fredholms theorem to pure surface-nonlinearities? 
If yes, what is the adapted route for defining a linear adjoint operator. If not, what is a generic and mathematically 
proven alternative since the hitherto used routes are lacking these features. Therefore further analyses need to be 
done in order to solve this problem and to clarify the open questions. 

If these problems for the Rosensweig instability, caused by the single external excitation of a magnetic field, are 
solved, other instabilities caused by different external excitations can be fruitful tackled. Phenomena as standing 



twin peaks [|39| or domain structures | 
designated future examples for an ana 



J] for parametrically excited MF under the influence of a magnetic field are 
ysis by amplitude equations derived from the basic hydrodynamic equations. 
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